The Range of Validity for the Kelvin Force 
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In a recent Letter Luo, Du and Huang reported a 
novel convective instability driven by a force rarely stud- 
ied before - that exerted by an external magnetic field on 
a strongly magnetizable liquid. The associated physics is 
surprisingly rich and promises many more interesting re- 
sults for the future. Unfortunately, the analysis starts 
from a misconception and employs the Kelvin force out- 
side its range of validity. Since few would recognize this 
as a mistake, and since its consequence in the given ex- 
periment is particularly direct and critical, this is a point 
well worth being clarified, and clearly understood. 

In the experiment, ferrofluid is exposed to a constant 
B-field. Yet, since the temperature T and the density 
p of magnetic particles vary, so does the magnetic field 
H = B/[l + x(T, p)], giving rise to a finite Kelvin force. 
With x the magnetic susceptibility and M = xH the 
magnetization, this force is given as 


f = MiVHi = (M t B t )V 
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(1) 


(Summation over the index i is implied.) Eq (1) may be 
derived from the more general Helmholtz force || , 


-V{\H 2 pdx/dp) - \H 2 V X , 


(2) 


by considering a dilute ferrofluid, and taking x as pro- 
portional to the particle density p, or pdx/dp — x- Then 
Eq (2) clearly reduces to f = jxV{H 2 ) = MiVHi. 

All this seems rather convincing, but in fact hides a 
pitfall. Closer scrutiny reveals that f = MiVHi is only 
valid to linear order in \. (Except in unconventional sys- 
tems of more recent dates, the magnetic susceptibility \ 
is usually much smaller than 1 , so terms of higher order 
in x have always been negligible. This may well be the 
reason why the confined range of validity of the Kelvin 
force has been such a well kept secret.) If true, the ex- 
pressions of Eq (1) merely states that the force vanishes 
- to linear order in \. No result derived from Eq (1) is 
then trustworthy. 

To qualitatively understand this restriction, define a 
different susceptibility, M = xB- With the permeability 
given as^ = l + x = (l- x) _ \ w e have x = x/(l + x)- 
Both susceptibilities are clearly physically equivalent, 
and we have no a priori reason to prefer either. Em- 
ploying dx = dx/(l + x) 2 : we may rewrite Eq (||) as 

f = +V(\B 2 Pa dx/d Pa ) - \B 2 V X - (3) 

This time, assuming x as proportional to p, we obtain 


f = MiVBi, 


(4) 


a result obviously different from Eq (1) - but one that 
also vanishes for uniform B-fields, so there is no disagree- 
ment to f = MiVHi in linear order. 


Now, since Eqs (||) and (||) are algebraically equiva- 
lent, the difference must lie between the two seemingly 
innocuous assumptions, x or x ~ P- Reviewing the above 
derivations, it is obvious that if one of the two assump- 
tions were strictly correct, the other would be wrong, and 
only the associated force expression is applicable. 

Generically, on the other hand, both x and \ are power 
series of p. So we are simply approximating, discarding 
quadratic and higher order terms, when we assume that 
either is linear in p. The consistent dilute limit is given 
when all terms ~ p 2 (and higher) are discarded. With 
X ~ p, this necessarily implies that we must also discard 
all terms ~ x 2 - As a result, x = x/(l + x) ~ X and 
MiVBi w MiVHi. We conclude: The Kelvin force is 
valid to linear order in the density p and the susceptibility 
X (or magnetization Mi). Especially, both MiVBi and 
MiVHi are valid expressions for the Kelvin force. 

The force in the experiment of |jj is of course finite. 
A proper, quantitative evaluation is given by including 
terms of higher order in p. As a first-step, we consider 
the next order terms: 


X 


X = ap(l + f3ap-\ ), 

ap[l + {13- l)ap-\ ]. 


(5) 
(6) 


(Note that with XiX = a P m the dilute limit, ap is 
simply the sum of single particle contributions, from 
non-interacting dipoles.) Inserting these expansions into 
Eq (|^) or (g), we find, for a constant £L field, 


f = ±B 2 V[(/3-l)(ap) 2 


(7) 


There are four different microscopic models usually 
employed to calculate the temperature and concentration 
dependence of the susceptibility of ferrofluids: Weiss, On- 
sager, mean-spherical, and high-temperature. Though 
different in details, all provide the same value (3 = 1/3, 
in good agreement with experimental data pi. 

Assuming x ~ p or x ~ p to hold strictly is respectively 
equivalent to f3 = and [3 = 1, with additional restric- 
tions for the yet higher order terms. Both assumptions 
are arbitrary, and in stark contrast to our microscopic 
understanding of magnetism. 
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